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Abstract 

The resonance states of one- and two-particle Hamiltonians are studied using variational expan- 
sions with real basis-set functions. The resonance energies, E^, and widths, F, are calculated using 
the density of states and an C"^ golden rule-like formula. We present a recipe to select adequately 
some solutions of the variational problem. The set of approximate energies obtained shows a very 
regular behaviour with the basis-set size, N . Indeed, these particular variational eigenvalues show 
a quite simple scaling behaviour and convergence when N ^ oo. Following the same prescription 
to choose particular solutions of the variational problem we obtain a set of approximate widths. 
Using the scaling function that characterizes the behaviour of the approximate energies as a guide, 
it is possible to find a very good approximation to the actual value of the resonance width. 
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I. INTRODUCTION 



The methods used to calculate the energy and the lifetime of a resonance state are nu- 



merous IllQ and, in some cases, has been put forward over strong foundations 



Bl. However, 



the analysis of the numerical results of a particular method when applied to a given problem 
is far from direct. The complex scaling (complex dilatation) method when the Hamil- 
tonian H allows its use, reveals a resonance state by the appearance of an isolated complex 
eigenvalue on the spectrum of the non-Hermitian complex scaled Hamiltonian, H{6) jlj]. Of 
course in an actual implementation the rotation angle 9 must be large enough to rotate 
the continuum part of the spectrum beyond the resonance's complex eigenvalue. Moreover, 
since most calculations are performed using finite variational expansions it is necessary to 
study the numerical data to decide which result is the most accurate. To worsen things 
the variational basis sets usually depend on one (or more) non-linear parameter. For bound 
states the non-linear parameter is chosen in order to obtain the lowest variational eigenvalue. 
For resonance states things are not so simple since they are embedded in the continuum. 
The complex virial theorem together with some graphical methods [o], [l^ allows to pick the 



best numerical so^ 
the 6 trajectories 



ution o 
9 



101- 



a given problem, which corresponds to the stabilized points in 



Other methods to calculate the energy and lifetime of the resonance, based on the numer- 
ical solution of complex Hamiltonians, also have to deal with the problem of which solutions 
(complex eigenvalues) are physically acceptable. For example, the popular complex absorb- 
ing potential method, which in many cases is easier to implement than the complex scaling 
method, produces the appearance of nonphysical complex en ergy stabilized points that must 



be removed in order to obtain only the physical resonances 

The aforementioned issues explain, at some extent, why the methods based only in the 
use of real variational functions are often preferred to analyze resonance states. These 



techniques reduce the problem to the calculation of eigenvalues of real symmetric matrices 



a 



12 



Iq . Of course, these methods also have its own drawbacks. One of the main problems 



was recognized very early on (see, for example, the work by Hol0ien and Midtdal 16|): if the 
energy of an autoionizing state is obtained as an eigenvalue of a finite Hamiltonian matrix, 
which are the convergence properties of these eigenvalues that lie in the continuum when the 
size of the Hamiltonian matrix changes? But in order to obtain resonance-state energies it 
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is possible to focus the analysis in a global property of the variational spectrum: the density 
of states (D0S)[l3], being unnecessary to answer this question. 

The availability of the DOS allows to obtain the energy and lifetime of the resonance in a 
simple way, both quantities are obtained as least square fitting parameters, see for example 



6|, llJI- Despite its simplicity, the determination of the resonance's energy and width based 
in the DOS is far from complete. There is no a single procedure to asses both, the accuracy 
of the numerical findings and its convergence properties, or which values to pick between 
the several "candidates" that the method offers 

Recently, Pont ei a/ jl3| have used finite size scaling arguments 18 1 to analyze the prop- 
erties of the DOS when the size of the Hamiltonian changes. They presented numerical 
evidence about the critical behavior of the density of states in the region where a given 
Hamiltonian has resonances. The critical behavior was signaled by a strong dependence 
of some features of the density of states with the basis-set size used to calculate it. The 
resonance energy and lifetime were obtained using the scaling properties of the density of 
states. However, the feasibility of the method to calculate the resonance lifetime laid on the 
availability of a known value of the lifetime, making the whole method dependent on results 
not provided by itself. 

The DOS method relies on the possibility to calculate the Ritz-variational eigenfunctions 
and eigenvalues for many different values of the non-linear parameter rj (see Kar and Ho 
For each basis-set size, A^, used, there are variational eigenvalues E^ij])] m = 1, . . . , N . 
Each one of these eigenvalues can be used, at least in principle, to compute a DOS, p'^{E), 
resulting, each one of these DOS in an approximate value for the energy, Er ~ E^, and 
width, r ~ r^V' the resonance state of the problem. If the variational problem is solved 
for many different basis-set sizes, there is not a clear cut criterion to pick the "better" result 
from the plethora of possible values obtained. This issue will be addressed in Section lllli 

In this work, in order to obtain resonance energies and lifetimes, we calculate all the 
eigenvalues for different basis-set sizes, and present a recipe to select adequately certain 
values of A^, and one eigenvalue for each N elected, that is, we get a series of variational 
eigenvalues {E^^{r])}. 

The recipe is based on some properties of the variational spectrum which are discussed in 
Section [TTl The properties seem to be fairly general, making the implementation of the recipe 
feasible for problems with several particles. Actually, because we use scaling properties for 
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large values of A^, the applicability of the method for systems with more than three particles 
could be restricted because the difficulties to handle large basis sets. 

The set of approximate resonance energies, obtained from the density of states of a 
series of eigenvalues selected following the recipe, shows a very regular behaviour with the 
basis set size. This regular behaviour facilitates the use of finite size scaling arguments to 
analyze the results obtained, in particular the extrapolation of the data when — > oo. 
The extrapolated values are the most accurate approximation for the parameters of the 
resonance state that we obtain with our method. This is the subject of Section llVt where 
we present results for models of one and two particles. 

Following the same prescription to choose particular solutions of the variational prob- 
lem we obtain a set of approximate widths in Section |Vl Using the scaling function that 
characterizes the behaviour of the approximate energies as a guide, it is possible to find a 
very good approximation to the resonance width since, again, the data generated using our 
prescription seems to converge when N oo. Finally, in Section I VI I we summarize and 
discuss our results. 

II. SOME PROPERTIES OF THE VARIATIONAL SPECTRUM 

When one is dealing with the variational spectrum in the continuum region, some of its 
properties are not exploited to obtain more information about the presence of resonances, 
usually the focus of interest is the stabilization of the individual eigenvalues. The stabiliza- 
tion is achieved varying some non-linear variational parameter. If rj is the inverse charac- 
teristic decaying length of the variational basis functions, then the spectrum of the kinetic 
energy scales as 77^, moreover, for potentials that decay fast enough, the spectrum of the 
whole Hamiltonian also scales as 77^ for large (or small) enough values of 77 (see appendix). 
This is so, since the variational eigenfunctions are approximations to plane waves except 
when ?7 belongs to the stabilization region. When rj belong to the stabilization region of 
a given variational eigenvalue, say i?]^(?7), then E^{ri) ~ (where is the resonance 
energy) and the variational eigenfunction ip^irj) has the localization length of the potential 
well. We intend to take advantage of the changes of the spectrum when 77 goes from small 
to large enough values. 

The variational spectrum satisfies the Hylleras-Undheim theorem or variational theorem: 
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if is the basis set size, and is the m — th eigenvalue obtained with a variational basis 
set of size A^, then 

E^{r,)>E^^M- 

Actually, since the threshold of the continuum is an accumulation point, then for small 
enough values of t] and a given j G N there is always a /c G N such that 

small ) ■ (1) 

For the kinetic energy variational eigenvalues, and for fixed A^, m,j, and k, if the ordering 
given by Equation ([T]) holds for some t] then it is true for all 77. Of course this is not true 
for a Hamiltonian with a non zero potential that support resonance states. So, we will take 
advantage of the variational eigenvalues such that for 77 small enough satisfies Equation ([T]) 
but, for 7] large enough 

E'^iVlarge) < E'^'l{{'qiarge) ■ (2) 

Despite its simplicity, the arguments above give a complete prescription to pick a set of 
eigenstates that are particularly affected by the presence of a resonance. Choose A^ and m 
arbitrary, and then look for the smaller values of j and k such that the two inequalities. 
Equations. ([T]) and ([2]) are fulfilled. So far, all the examples analyzed by us show that if 
the inequalities are satisfied for some j and k then they are satisfied too by the eigenvalues 
^K^,forn = 2,3,4,.... 

III. MODELS AND METHODS 

To illustrate how our prescription works we used two different model Hamiltonians. The 



first model, due to Hellmann 



2l|, is a one particle Hamiltonian that models a A"e-electron 



ow energy 



13|. 



atom. The second one is a two particle model that has been used to study the 
and resonance states of two electrons confined in a semiconductor quantum dot 

The details of the variational treatment of both models will be kept as concise as possible. 
The one particle model has been used before for the determination of critical nuclear charges 
for A"e-electron atoms [22 1, it also gives reasonable results for resonance states in atomic 



anions 



19| and continuum states 23|. The interaction of a valence electron with the atomic 



core is modeled by a one-particle potential with two asymptotic behaviours. The potential 
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behaves correctly in the regions where electron is far from the atomic core (A'e — 1 electrons 
and the nucleus of charge Z) and when it is near the nucleus. The Hamiltonian, in scaled 
coordinates r — ?■ Zr, is 

H = -iv^ - i + 7 (1 - ^-*) . (3) 
where 7 = { '^"z^ ) 5 is a range parameter that determines the transition between the 
asymptotic regimes, for distances near the nucleus Vo(r) ^ — ^ and in the case r — )■ 00 the 
nucleus charge is screened by the A'^e — 1 localized electrons and Voo{r) ^ ~ ^~zr~^^ ■ 

Another advantage of the potential comes from its analytical properties. In particular 
this potential is well behaved and the energy of the resonance states can be calculated using 
complex scaling methods. So, besides its simplicity, the model potential allows us to obtain 
the energy of the resonance by two independent methods and check our results. 

The two particle model that we considered describes two electrons interacting via the 
Coulomb repulsion and confined by an external potential with spherical symmetry. We use 
a short-range potential suitable to apply the complex scaling method. The Hamiltonian H 
for the system is given by 

H = -^Vl^ - ^Vl^ + V{n) + V{r,) + (4) 
2m ^ 2m ^ — ri| 

where Vir) = — (Vo/rQ) exp(— r/rg), Tj the position operator of electron i = 1,2; Tq and 
Vq determine the range and depth of the dot potential. After re-scaling with tq, in atomic 
units, the Hamiltonian of Equation (j4j) can be written as 

H = - I'^l - ^oe- - K,e- + (5) 



where \ = vq. 



le necessary 



The variational spectrum of the two particle model. Equation ([5]), and all t 
algebraic details to obtain it, has been studied with great detail in Reference 13[ so, until 
the end of this Section, we discuss the variational solution of the one particle model given 
by Equation ([3]). 

The discrete spectrum and the resonance states of the model given by Equation ([3]) 
can be obtained approximately using a real truncated basis set {^i{r])}^ to construct a 
NxN Hamiltonian matrix Hijijf) = ($4(77) | H |$j(?7)). We use the Rayleigh-Ritz Variational 
method to obtain the approximations {"ipNiv)) 
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N 



IV'M) = 5^cf)(r/)|<|.,(r/)) , m = l,---,N. (6) 



i=l 



For bound states this functions are variationally optimal. The functions |$j(?7)) are 

3/2 -r,r/2 

\HV)) = ^^——Lfivr) (7) 

V(2 + + 2) 

and Lf\rir) are the associated Laguerre polynomials of 2^*^ order and deCTee i. The non- 
linear parameter rj is used for eigenvalue stabilization in resonance analysis [71, ll5| • Note 
that T] plays a similar role that the finite size of the box in spherical box stabilization 
procedures 24], as stated by Kar et. al. 



Resonance states are characterized by isolated complex eigenvalues, E^-es = Er—iT/2, T > 
0, whose eigenf unctions are not square-integrable. These states are considered as quasi- 
bound states of energy Er and inverse lifetime F. For^he Hamiltonian Equation Q, the 
resonance energies belong to the positive energy range 



Using the approximate solutions of Hamiltonian ([3]) we analyze the DOS method [17 1 
that has been used extensively to calculate the energy and lifetime of resonance states, in 
particular we intend to show that 1) the DOS method provides a host of approximate values 
whose accuracy is hard to assess, and 2) if the DOS method is supplemented by a new 
optimization rule, it results in a convergent series of approximate values for the energy and 
lifetime of resonance states. 



The DOS method 



The DOS method relies on the possibility to calculate the Ritz- variational eigenfunctions 
and eigenvalues for many different values of the non-linear parameter rj (see Kar and Ho 
The localized DOS p{E) can be expressed as [l7| 



dE{r];-f,6) 



dr] 



(8) 



Since we are dealing with a numerical approximation, we calculate the energies in a dis- 
cretization {r]i} of the continuous parameter rj. In this approximation. Equation ([8]) can be 
written as 



p^{E) = p{E^{rj,-j,6)) 



E^\r]i^,;^,6)-E^{r],^,;-f,6) 



(9) 
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where E'^{rj; 7, 5) is the j'-th eigenvalue of the N x N matrix Hamiltonian with 7 and 5 
fixed. 



In complex scaling methods the Hamiltonian is dilated by a complex factor a 



As was pointed out long time ago by Moiseyev and coworkers 
and |a| are equivalent, in fact, our parameter rj corresponds to a 
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\a\ e 



-id 



the role played by 77 
0). Besides, the DOS 



attains its maximum at optimal values of rj and Er that could be obtained with a self- 
adjoint Hamiltonian without using complex scaling methods 9|. So, locating the position 
of the resonance using the maximum of the DOS is equivalent to the stabilization criterion 
used in complex dilation methods that requires the approximate fulfillment of the complex 
virial theorem 10|. 



The values of Er{'~f,S) and r{'~f,6) are obtained performing a nonlinear fitting of p{E), 
with a Lorentzian function, 



PiE) = po 



A 



r/2 



(10) 



7r[(E-E,)2 + (r/2)2]- 
One of the drawbacks of this method results evident: for each pair 7, 6 there are several 

p^{E), and since each p^{E) provides a value for Er{N, m) and r(A^, m) one has to choose 
which one is the best. Kar and Ho 3] solve this problem fitting all the p^{E) and keeping 
as the best values for Er and F the fitting parameters with the smaller value. At least 
for their data the best fitting (the smaller x^) usually corresponds to the larger j. This fact 
has a clear interpretation, if the numerical method approximates E^ with some Er{N,j), 
where is the basis set size of the variational method, a large j means that the numerical 
method is able to provide a large number of approximate levels, and so the continuum of 
positive-energy states is "better" approximated. 

In a previous work 13| we have shown that a very good approximation to the energy of 
the resonance state is obtained considering just the energy value where p^(-E) attains its 
maximum. We denote this value as Ej.{N,m). Figure [1] shows the approximate resonance 
energy Er{N, m) for different basis set size A^, where m is the index of the variational eigen- 
value used to calculate the DOS. We used the values 7 = 1.125 and 5 = 0.211 corresponding 
to the ones used before 19[ in the analysis of O resonances. The Figure [1] also shows the 
value calculated using complex scaling. It is clear that the accuracy of all the values shown 
is rather good (all the values shown differ in less than 6xl0~^), and that larger values on 
n provide better values for the resonance energy. These facts are well known from previous 
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FIG. 1: (color on-line) Er{N,m) vs 1/N for (7 = 1.125,(5 = 0.211). The values were obtained from 
the peak of the DOS resulting from Equation ([9]) for several basis set sizes, N = Nq + 3A;; k = 
0, . . . , km- The (orange) dots data correspond to A'^o = 50 m = 17 and km = 33; the (green) square 
dots to Nq = 51, m = 10 and km = 32; the (cyan) diamond dots to Nq = 51 m = 5 and km = 17; 
and the (violet) triangle dots to Nq = 51 m = 3 and km = 17. The cross shows the value we 
obtained using complex scaling {6 = it/5, N = 70). 

works, i.e. almost all methods to calculate the energy of the resonance give rather stable 
and accurate results for Er- However, the practical importance of this fact is reduced: these 
are uncontrolled methods, so the accuracy of the values obtained from the DOS can not be 
assessed (without a value independently obtained) and these values do not seem to converge 
to the value obtained using complex scaling when is increased and n is kept fixed. 

There is another fact that potentially could render the whole method useless: for small 
or even moderate m, the values Er{N,m) become unstable (see Figured]) when is large 
enough. This last point has been pointed previously [l^. In the problem that we are 
considering is rather easy to obtain a large number of variational eigenvalues in the interval 



where the resonances are located, allowing us to calculate K 



situation is far from common see, for example. References 



iV, m) up to m = 17, but this 



12 



13 



26|. 
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IV. SCALING OF THE RESONANCE ENERGY 



So far we have presented only results about the behaviour of the one particle Hamiltonian, 
from now on we will discuss both models, Equations and El 




FIG. 2: (color on-line) (a) Variational spectrum obtained for several basis set sizes of the one- 
particle Hamiltonian. There are two noticeable features associated to changes in the the density 
of curves. One of these features corresponds to the stabilization zone, where the derivative of 
the eigenvalues is minimal and, correspondingly, the density of states is larger. The other one 
corresponds to the gaps in the density of eigenvalues when E is kept fixed, (b) Detailed view of 
the data shown in panel (a). There is a well defined crossing between several eigenvalues, which 
form a bundle of states. The bundle is formed by £^5^, £^5!, ^gQ, -Eg|, and SgJ. 

It is known that the variational eigenvalues E^{ri) do not present crossings when they 
are calculated for some fixed values of N, i. e the variational spectrum is non-degenerate for 
any finite Hamiltonian matrix. As a matter of fact the avoided crossings between successive 
eigenvalues in the variational spectrum are the watermark of a resonance. An interesting 
feature emerges when the variational spectrum for many different basis set sizes N are 
plotted together versus the parameter 77. Besides the places where — attains its minimum 
value, which correspond to the stabilization points, there are some gaps which correspond to 
crossings between eigenvalues obtained with different basis set sizes, see Figure [2J Moreover, 
the crossings corresponds to eigenvalues with different index m, and are the states that satisfy 
the inequalities Equations ([T]),and ([2]). 
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It is worth to remark that the main features shown by Figure [2] are independent of the 
number of particles of the Hamiltonian and the particular values of the threshold of the 
continuum. Figure E] shows the behaviour of the variational eigenvalues obtained for the 
two particle Hamiltonian Equation In this case the ionization threshold is not the 

asymptotic value of the potential, but it is given by the energy of the one particle ground 
state. The resonance state came from the two-particle ground state that becomes unstable 
and enters into the continuum of states when the quantum dots becomes "too small" to 
accommodate two electrons. For more details about the model, see reference 13|. 

The left panel of Figure S] shows the behaviour of the maximum value of the DOS, 
Pmaxi^i"^)-! o'^^ particle Hamiltonian, obtained for different basis-set sizes and fixed 

m (in this case m = 17), and the pmax{N,m) obtained choosing a "bundle" of states that 
are linked by a crossing, these states have = 51, 54, 57, . . . , 168 and m = 16, 17, 18, . . . , 55 
respectively. From our numerical data, the maximum value of the DOS scales with the basis- 
set size following two different prescriptions. For m fixed, Pmax{N,m) ~ A^", with a > 0, 
while when the pair (A^, m) is chosen from the set of pairs that label a bundle of states 
Pmax{N,m) ~ A^"'^, with /3 > 0. In particular, for m = 17 we get that Pmax{N, 17) ~ j\r0.759^ 
and Pmax{N, m) ~ A^-^ when (A^, m) = (51, 17), (54, 18), (57, 19), . . .. 

Of course we can pick sets of states that are not related by a crossing. For instance, we 
also picked sets with a simple prescription as follows: choose a given initial pair (A'^o,mo) 
and form a set of states with the states labeled by (A'^o, mo), (A'^o + AA^, mo + Am), (Nq + 
2AN,mo + 2Am) and so on. Figure IH^a) shows two examples obtained choosing A'^o = 51, 
AA^ = 3 and mo = 16 and mo = 18 both with Am = 1. Quite interestingly, the data in 
Figure m show that the scaled maxima of the DOS for a bundle and two different sets seem 
to converge to the same value when A^ — oo, but only for the bundle the scaling function is 
A^~^. The advantage obtained from picking those eigenvalues E^^t]) that belong to a given 
bundle is still more evident when the corresponding DOS and Er are calculated. The right 
panel of Figure H] shows the Er obtained from the DOS whose maxima are shown in the left 
panel. It is rather evident that these values now seem to converge, besides, the extrapolation 
to A^ — )■ oo results in a more accurate approximate value for Er- In contradistinction, the 
values for E^ corresponding to a fixed index m (the values shown in the Figure H] correspond 
to m = 17) do not seem to converge anywhere close to the value obtained using complex 
rotation. 
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FIG. 3: (color on-line) Rayleigh-Ritz energies as a function of the nonlinear parameter rj. The 
calculations are for the resonances of the two-electron quantum dot model used in [l^. The 
Energies are calculated for Vq = 3 and A = 1.1 of equation ([5]). The dot and squares denote two 
different bundles both calculated with basis-set sizes Nq = 8; AA^ = 2. The eigenvalue numbers 
for the dotted bundle are given by rriQ = 4; Am = 1. For the bundle with squares the values are 
niQ = 3; Am = 1. The eigenvalues corresponding to the minimum (N = 8) and maximum (N = 20) 
basis-set size are single out. The dash-dotted cyan line is the resonance energy calculated using 
complex scaling. The dashed red lines indicates the exact energies where the continuum begins. 

Figure |5] shows the resonance energies obtained from the bundles of states shown in 
Figure [3] for the two-particle model. Since the numerical solution of this model is more 
complicated than the solution of the one-particle model the number of approximate values 
is rather reduced. However, it seems that the data also supports a linear scaling of Er{N, m) 
with 1 /N. 
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FIG. 4: . (a) Scaling of the DOS peak vs. the inverse basis-set size A^. One particle Hamiltonian. 
The figure shows the peak values with increasing N = 51, 54, 57, . . . , 168 , for m = 17 fixed (red 
crosses). The peaks for a bundle of eigenvalues (A'^o = 51; AA^ = 3; niQ = 17; Am = 1) show a 
linear scaling (black dots). For series of states that form bundles not associated to a crossing the 
scaling with 1/A^ is greater ( Nq = 51; AA^ = 3; mo = 16; Am = 1, green triangle up dots) or 
lower (A'o = 51; AA^ = 3 ; niQ = 18 ; Am = 1, orange triangle down dots) than 1. (b) The energy 
position of the DOS peaks as a function of 1/A^. The figure shows that the fixed m = 17 resonance 
energies (red crosses) do not converge to the resonance energy calculated by complex scaling (bold 
black cross). The scaling for the bundle, and near bundle, series of states is almost linear and an 
extrapolated value can be estimated. 

V. SCALING OF THE RESONANCE WIDTH 

Many real algebra methods to calculate resonance energies use a Golden Rule-like formula 
to calculate the resonance width. In this section we will use the formula and stabilization 
procedure proposed by Tucker and Truhlar jsol that we will describe briefly for completeness. 



This projection formula seems to wor 



i better for one-particle models. For two-particle 



models its utility has been questioned 31|, so to analyze the width of the resonance states 



of the quantum dot model we fitted the corresponding DOS using Equation (11 01) . 



The method of Tucker and Truhlar 30|] is implemented by the following steps. Choose 
a basis {(f) jif])} where ?7 is a non-linear parameter. Diagonalize the Hamiltonian using up 
to functions of the basis. Look for the stabihzation value Sres and its corresponding 
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FIG. 5: Resonance Energy for the two bundles (dots and squares) showed in Figure [3l Since there 
is an even-odd behavior the points of each bundle are grouped in two sets. Data in sets 1 ( Black 
) correspond to = 8, 12, 16, 20, and for sets 2 (green) correspond to = 10, 14, 18. The black 
cross shows the width calculated using complex scaling. 

eigenfunction ipres which are founded for some value r]^. Define the projector 

= 1^-'^) (^-.''1 , (11) 

where \ip^'^) is the normalized projection of tpres onto the basis {(pj{v)}j'=i for other 77. 

Diagonalize the Hamiltonian H = {1 — Q^)H{1 — Q"^) in the basis {(pj{v)}j'=i^ again as a 
function of r], and find a value t]'^ of r] such that 

Enivl = Sres, (12) 

where En{ri'^) denotes eigenvalue n of the projected Hamiltonian for the scale factor t]'^, and 
Xniv'^) is the corresponding eigenfunction. 

With the previous definitions and quantities, the resonance width F is given by 

F = 27rp(r/^£,es) \{Aes \H\ Xn{vl)f , (13) 

where 

p(r/^£,e.) = 2/ [En+iir]") - En^,{r]')] . (14) 
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FIG. 6: Scaling of the width vs. the inverse basis-set size N. All points correspond to the same 
bundle of states {N = 51, 54, 57, ... , 150 ; m = 17, 18, 19, ... , 50). The curves were obtained with 
different values of the projection eigenvalue n ( 1, 2, 3, 4, 5, from bottom to top) . The black line is 
a fit for n = 3 (Equation (jl7p ) and the lowest n curve that shows a scaling exponent greater than 
one (7(3) ~ 1.2). 

Despite some useful insights, the procedure sketched above does not determine all the 
intervening quantities, for instance there are many solutions to equation f|T2|) and, of course, 
the stabilization method provides several good candidates for tpres and r]^. 

We are able to avoid some of the indeterminacies associated to the Tucker and Truhlar 
procedure using a bundle of states associated to a crossing, so tpres and r]'' are given by 
any of the eigenf unctions associated to a bundle and ri^'{N^m) comes from the stabilization 
procedure. Then we construct projectors 

Ql,m=\n{ri)){n{r,)\, (15) 

where il''p}{v) is one of the variational eigenstates that belong to a bundle of states. With 
the projectors QJ^^ construct Hamiltonians H, and find the solutions to the problem 

Eniv:) = Sres. (16) 

Since there is not an a priori criteria to choose one particular solution of Equation (fT6|) we 
show our numerical findings for several values of n. Figure |6] shows the behaviour of the 
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resonance width calculated with Equation ( !T3|) . where we have used ipll, i^lf, i^lj, • • • as ipres 
and Xn, where n = 1, 2, 3, 4 and 5. 

Despite that the different sets corresponding to different values of n do not converge to 
any definite value, for N large enough all the sets scale as N~^^"'\ with 7(n) > 0. Since the 
resonance energy scales as A^~^, at least when a bundle of states with a crossing is chosen 
to calculate approximations (see Figure |4]), we suggest that the right scahng for T is given 
by 7 = 1. Of course for a given basis size, particular variational functions, stabilization 
procedures and so on, we can hardly expect to find a proper set of F whose scaling law 
would be A^~^. Instead of this we propose that the data in the right panel can be fitted by 

F;v,n ~ A{-f{n)) + 5(7(r2))iV-TW, ^{n) > 0. (17) 

Then the best approximation for the resonance width is obtained fitting the curve and 
selecting the F^ as A{'y{n)) for j{n) the closest value to one. 
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FIG. 7: (color on-line) Resonance widths for the two bundles (dots and squares) showed in Figure[3l 
Since there is an even-odd behavior the points of each bundle are grouped in two sets. As in figure 
[5l data in sets 1 ( Black ) correspond to = 8, 12, 16, 20, and for sets 2 (green) correspond to 
N = 10,14,18. The widths where calculated using the method described in Reference 71] . The 
black cross is the width calculated using complex scaling. 



As pointed in Reference 



31| , the projection technique to calculate the width of a resonance 
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can be implemented if a suitable form of the projection operator can be found. As this 
procedure is marred by several issues we used the DOS method to obtain the approximate 
widths of a resonance state of the two particle model. Figure [7] shows the widths calculated 
associated to the energies shown in Figure the parameters of the Hamiltonian are exactly 
the same. 

There is no obvious scaling function that allows the extrapolation of the data but, even 
for moderate values of A^, it seems as the data converge to the value obtained using complex 
scaling. 

VI. CONCLUSIONS 

In this work we analyzed the convergence properties of real C? basis-set methods to obtain 
resonance energies and lifetimes. The convergence of the energy with the basis-set size for 
bound states is well understood, the larger the basis set the better the results and these 
methods converge to the exact values for the basis-set size going to infinite (complete basis 
set). This idea is frequently applied to resonance states. The increase of the basis-set size in 
some commonly used methods does not improve the accuracy of the value obtained for the 
resonance energy , as showed in Figure [H This undesirable behavior comes from the fact 
that the procedure is not variational as in the case of bound states. Moreover, the exact 
resonance eigenfunction does not belong to the Hilbert space expanded by the complete basis 
set. In this work we presented a prescription to pick a set or bundle of states that has linear 
convergence properties for small width resonances. This procedure is robust because the 
choice of different bundles results in very similar convergence curves and energy values. In 
fact, in the method described here, the pairs (A^, m) of the bundles play the role of a second 
stabilization parameter together with the variational parameter 77. of a second We tested 
the method in others one and two particle systems and the general behavior of them is the 
same. The results are very good in all cases leading to an improvement in the calculation of 
the resonance energies. Nevertheless we have to note that the method could no be applied in 
cases where two or more resonance energies lie very close because the overlapping bundles. 

The lifetime calculation is more subtle. The use of golden-rule-like formulas, as we applied 
here, always give several possible outcomes for the width FAr^n; corresponding to different 
pseudo-continuum states |xn)- The projection technique. Equation (IT^ . is not the exception 
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and it is not possible to select a priori which value of {Fat is the most accurate. The linear 
convergence of the DOS with basis-set size suggests that the scaling in the lifetime value, in 
accordance with the Energy scaling, should be linear. Regrettably, the projection method 
gives discrete sets of values which cannot be tuned to obtain an exact linear convergence. 
Our recipe is to choose the set r7v,n. whose scaling is closest to the linear one, then the best 
estimation for the resonance width is obtained from extrapolation. 

Many open questions remain on the analysis of the different convergence properties of 
resonance energy and lifetime. The method presented here to obtain the resonance energy 
from convergence properties works very well, but the appearance of bundles in the spectrum 
is not completely understood. Even there is not a rigorous proof, the numerical evidence 
supports the idea that the behaviour of the systems studied here is quite general. 



Appendix 

In this appendix we give arguments that support our assumptions on the scaling of 
the eigenenergies with the basis-set parameter 77. We present our argument for one body 
Hamiltonians, but it is straightforward to generalize to more particles with pair interactions 
decaying fast enough at large distances. 

First et all, we present two very known results that we need later. 

1. Let A be an n X n matrix with all its matrix elements having the form aij{T]) = 
ajj(l) f{r]), where /(I) = 1, then if /(r/) 7^ the eigenvalues of A scales with f{r]): 

det[Aiv) - Xiv) /] = ^ det[Ail) - M j] = ^ A(l) = ^ . (A.18) 

2. Let A, B n X n he symmetric matrices with \bij\ < const., and Aj, i = 1, ■ ■ ■ , n 
the eigenvalues of A, B and A + eB respectively, in nondecreasing order, then, by the 
minimax principle 32 1 



ai + el3i<\i<a^ + e[3n (A.19) 

Consider a spherical one-particle potential with compact support, t;(r) = if r > i?, and 
finite, \v{r)\ < 00 Vr (both conditions could be relaxed, but we adopt them for simplicity). 
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Let the basis-set functions be of the form $„(?7,r) = c„(?7) ^niv^) with ||$„(?7, r)|| = 1 and 
Cn(l) = 1, then the coefficients take the form Cn{f]) = Cn{l)r]'^^^- 

With these assumptions, the matrix elements for the kinetic energy are 



(A.20) 



c*Jl)c^{l)v^ J d3a;$„(l,r)^|: (r^f) $„(1, r) = r^'T^A^) , 
and then, by Equation ( 1A.18I) . all the eigenvalues of the kinetic energy have the same scaling 
with rj"^. We have to show that, in both limits, 77 —t- and r/ — )■ 00, for all the potential 
matrix elements hold Vm,n/v'^ ^ 1; ^ind then, by the Wielandt-Hoffman theorem 32|, the 
eigenenergies are a perturbation of the eigenvalues of the kinetic energy. 



77 



In this limit 77/? < 1 ^{rjr) ~ $(0) for re(0, R], then 



47r<l>;;(l,0)<l'(l,0) / r^drv{r) 
Jo 



where Vm.,n do not depend on rj. Then, in this limit Vm,n/v'^ ~ r/ ^ 1. 



7] ^ 00 



For ?7 i? — )■ 00 we obtain for the potential matrix elements 



(A.21) 



(A.22) 



C(l)c„(l)^(0). 
Then also in this limit we obtain Vm,n/ff' ~ l/?]^ <^ 1. 
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